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Expansions of the Wilson determinant in lattice QCD with quarks produce gauge action terms which shift 
the coupling constant of the fundamental representation plaquette action and induce an adjoint representation 
plaquette action. We study the magnitude of these induced couplings with two flavors of Wilson fermions. We 
utilize a microcanonical demon method, which allows us to measure the induced couplings directly from gauge 
configurations generated by full fermionic simulations. 



1. Introduction 

Lattice QCD simulations with two light Kogut- 
Susskind quarks seem to indicate that QCD with 
two massless quarks has a second order finite tem- 
perature deconfining phase transition, which is 
transformed into a smooth crossover if the quarks 
are massive. This scenario is also consistent with 
the universality arguments [1]. However, with two 
fiavors of Wilson quarks recent simulations have 
shown that at large values of the hopping param- 
eter the phase diagram is more complicated [2]: 
with Ni = 6 and at hopping parameter k ~ 0.19 
the system has a strong first order like transi- 
tion, where the plaquette expectation value has 
a sharp discontinous jump. However, the aver- 
age Polyakov loop remains small at this point, 
and starts to increase only when /3 is consider- 
ably larger. This behaviour is strongly reminis- 
cent of separated bulk and thermal transitions; 
this occurs, for example, in pure gauge SU(3) sim- 
ulations with mixed fundamental - adjoint action 
[3]. In fact, these systems resemble each other so 
much that one is tempted to assume that the dy- 
namical reason for these transitions is the same; 
i.e. that the Wilson fermionic action induces a 
strong adjoint pure gauge action. 

The pure gauge SU(3) fundamental-adjoint ac- 
tion can be written as S = PfSf -\- PaSa, where 
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where Tr Up is the fundamental representation 
trace. The gauge action induced by the fermion 
determinant can be studied perturbatively via the 
hopping parameter expansion; the leading terms 
are A/3p ex «■* and A/3^ ex k^'^ . This expansion is 
accurate only when k is very small. Results with 
much broader validity range can be obtained with 
the heavy quark perturbation. This method was 
used by Hasenfratz and DeGrand [7] to calculate 
A/3p, which agrees very well with the MC data. 
However, this expansion has not been performed 
for A/3^. 

2. Demon Algorithm 

We use the microcanonical demon method [4- 
6] to project out the induced gauge action from 
full fermionic simulations. Our ansatz for the ef- 
fective action is given by eqs. (1,2), with a pri- 
ori unknown coupling constants (3i = (3f or (3a ■ 
For both of the coupling constants we introduce 
a demon , which is a real-valued action variable 
with < A < Df^^^. Effectively, the fermionic 
action is used as a heat bath to thermalize the 
demons: a configuration generated with the orig- 
inal action is updated microcanonically with the 
demons keeping Si + Di constant for all i sepa- 
rately. At the end of the microcanonical update, 
we discard the old configuration and substitute 
it with a new one, while preserving the values 
of the demons. Repeating this many times the 
demons attain an equilibrium distribution and we 
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Figure 1. The induced /Sp. Small crosses on 
horizontal dashed lines indicate the simulation 
(j3p,K), opaque circles and black dots the mea- 
sured iSp-values using the effective fundamental- 
adjoint and fundamental only gauge action, re- 
spectively. The length of the vertical bars gives 
the magnitude of the induced A/Sp. 



Figure 2. The induced /3^. The scale on the right 
gives the magnitude of the induced i3a, measured 
as the vertical distance between the plot sym- 
bols and dashed horizontal lines. Points below 
the dashed lines indicate negative /3^. 



can measure (Di). The induced l3i can be solved 
from the equation 
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In our case, only the adjoint action demon was 
bounded from above (D^^^ = cxd). 

There are various non-equivalent ways to per- 
form the demon update. For each starting con- 
figuration, one can update microcanonically until 
the demons and the system are properly thermal- 
ized, or one can stop the microcanonical update 
after one - or even partial - update sweep. If 
the effective action is completely equivalent to the 
original action, these two methods give the same 
unique answer. However, in this case we have a 
dramatically truncated effective action, and these 
method need not be equivalent. Nevertheless, the 
differences are minor if the heat capacity of the 
demons is much smaller than that of the system, 
and in all our tests the possible differences were 
completely overwhelmed by the statistical noise. 



3. Simulations and Results 

We tested the demon method by applying it 
to pure gauge fundamental-adjoint simulation. 
With 4'*-lattice, we simulated the system at 
couplings (j3f,I3a)o = (3.6,1.8) and (4.0,2.0). 
The induced couplings were (3.594(6), 1.812(17)) 
and (4.017(27), 1.92(6)), respectively, compatible 
with the input values. The latter coupling pair is 
very close to the bulk transition line in (Pf,i3a)- 
plane. When we used an effective action con- 
sisting only of the S'f-part, /Sp was 4.688(9) and 
6.044(12); the latter value is very close to the ex- 
tension of the bulk transition line to the /3^ = 
-axis [3]. 

If the Wilson fermion action induces a strong 
adjoint coupling giving rise to bulk fundamental- 
adjoint transition, one should be able to observe 
the induced couplings already in small volumes. 
We performed simulations on 4^* lattices with 28 
different (j3p,K) pairs. In figs. 1 and 2 we show 
the measured jSp and Pa calculated with jSp = 
4.0, 4.9, 5.46, and several k values; in fig. 3 with 
constant k = 0.19. In figs. 1 and 3 we also show 
l3p when Seff = PpSp only. 
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Figure 3. The induced /Sp and i3a, when k = 0.19. 



When K is small, the quarks are very massive 
and the induced couplings are quite small (left 
side of figs. 1 and 2). When k is increased, we 
approach the critical line where niq -^ 0, and 
the fermionic contribution to the action becomes 
more significant. This is clearly visible as an in- 
crease in i3f in figs. 1 and 3. The critical values 
of K are approximately 0.16 (j3p = 5.46), 0.19 
(4.9) and 0.22 (4.0). We observe no significant 
increase in i3a ', on the contrary, when k^ is ap- 
proached, induced (3a becomes slightly negative! 
The minor role of the adjoint action is also ev- 
ident from the fact that (3f remains practically 
the same whether we use the 5^ -term of the ef- 
fective action or not. We also checked the results 
with a few simulations on 6^* lattices with similar 
results. 

In fig. 4 we compare A/3i? to the predictions 
of ref. [7], as a function of quark mass Wqa = 
K~^—K~^. The agreement is very good, especially 
when P% = 4.9. 

To conclude with, we did not observe that Wil- 
son fermions induce any significant adjoint pure 
gauge action, and that the induced fundamental 
action is very well described by analytical calcu- 
lations. It is very improbable that the transi- 
tion observed in #( = 6 Wilson thermodynam- 
ical simulations [2] could be explained by the 
fundamental-adjoint pure gauge transition, and 



Figure 4. Comparison of the analytical result for 
A/Sp [7] by Hasenfratz and DeGrand (solid line) 
to the MC data. 



the real cause of this phenomenon remains to be 
uncovered. 
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